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, ^ , Abstract 

I We consider the classical superintegrable Hamiltonian system given by 



> 

00 ^^ = '^ + ^=2(l + Aq2) + 2(l + Aq2)' 

m 

CO where W is known to be the "intrinsic" oscillator potential on the Darboux spaces of 

nonconstant curvature determined by the kinetic energy term T and parametrized by 
^-H A. We show that Hx is Stackel equivalent to the free Euclidean motion, a fact that 

directly provides a curved Pradkin tensor of constants of motion for H\. Furthermore, 
we analyze in terms of A the three different underlying manifolds whose geodesic motion 
J> is provided by T. As a consequence, we find that Hx comprises three different nonlinear 

physical models that, by constructing their radial effective potentials, are shown to be two 
different nonlinear oscillators and an infinite barrier potential. The quantization of these 



two oscillators and its connection with spherical confinement models is briefly discussed. 
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1 Introduction 



Let us consider the A^-dimensional (A^D) classical Hamiltonian defined by 

%=r(q,p)+Z.(q) = 5^^ + 5(^, (1) 

where A and u; are real parameters, and q, p G are conjugate coordinates and momenta 
with canonical Poisson bracket {qi,Pj} = Sij. 

The mathematical and physical relevance of this system rely on two main properties [1] : 
(i) Tix is a maximally superintegrable (MS) Hamiltonian, since it is endowed with the max- 
imum possible number of 2N — 1 functionally independent integrals of motion; and (ii) the 
central potential ^(q) can be interpreted as the "intrinsic" oscillator on the underlying curved 
manifold defined through the kinetic term T- In particular, T determines the geodesic mo- 
tion of a particle with unit mass on a conformally flat space which was constructed in [21 13] 
and is the A^D spherically symmetric generalization of the Darboux surface of type III [HIS]. 
The corresponding metric and scalar curvature depend on A and are given by 

^ 2 ^ , 2, , 2 , , (jV-l)(2jV + 3(iV-2)Aq2) 
ds =(1 + Aq)dq, i?(q) = -A (l + Aq^)^ 

From this viewpoint, Hx can be regarded as a MS "A-deformation" of the A^D isotropic 
harmonic oscillator with frequency uj since Huix^q'Hx = ^p^ + icu^q^. 

We recall that T-Lx can be identified as a particular case within other frameworks such as: 
(i) the "3D multifold Kepler" Hamiltonians [6, 7J (which generalize the MIC-Kepler and Taub- 
NUT systems); (ii) the "3D Bertrand systems" [U [9l [10] (coming from a generalization of 
the classical Bertrand's theorem [11] to curved spaces); and (iii) the "A^D position-dependent 
mass systems" [121 [131 ttH [151 [iHl [13 [El [13 [20] (see also references therein) provided that the 
conformal factor of the metric ^ is identified with the variable mass function m(q) = 1+Aq^. 

The aim of this paper is twofold. On one hand, in the next section we provide a deeper 
insight in the set of integrals of motion of T-Lx given in [1] by applying the so-called Stackel 
transform or coupling constant metamorphosis |21 1 1221 [23l (2^ [25j . In this way, we obtain the 
corresponding A-deformation of the Fradkin tensor of integrals of motion |26j for the isotropic 
harmonic oscillator. On the other hand, we explicitly show that T-Lx gives rise, in fact, to three 
different physical models. For this latter (and main) purpose, we present in section 3 which 
are the underlying manifolds that come out according to the values of A. This analysis leads 
to three types of manifolds which, in turn, correspond to two nonlinear oscillator systems plus 
a barrier-like one, which are studied in section 4 by constructing their associated effective 
potential. The final result is that the Hamiltonian T-Lx comprises the hyperbolic oscillator 
(A > 0), the spherical one (the "interior" space with A < 0) and an infinite potential barrier 
(the "exterior" space with A < 0). Remarkably enough, the effective oscillator potentials 
are, in this order, hydrogen-like and oscillator-like, which means that the quantization of Tix 
would provide different types of spherical confinement models like, for instance, [271128]. First 
results in this direction [29] are briefiy sketched. 
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2 Superintegrability and the Stackel transform 



The MS property of T-L\ is characterized by the fohowing statement. 

Theorem 1. (i) The Hamiltonian %\ for any real value of X, is endowed with the 
following constants of motion. 

• {2N — 3) angular momentum integrals: 

(j(m) ^ J2 {qiPj - qjPif , (^(m) = ^ {QiPj-QjPif, (3) 

l<«<i<m N-m<i<j<N 

where m = 2, . . . , N and C^^^ = Cf^jqy 

• A^^ integrals which form the ND curved Fradkin tensor: 

hj = PiPj - {2\'Hx{(l, P) - ^^)QiQj, (4) 
where i, j = 1, . . . , N and such that Tlx = ^ J2jLi ^a- 

(ii) Each of the three sets {Hx, C^}, {Hx, (m = 2, . . . , N ) and {I^^} (i = l,...,N) 

is formed by N functionally independent functions in involution. 

(Hi) The set {T-Lx}C^"^\C(^m)i hi} for m = 2,...,N with a fixed index i is constituted by 
2N — 1 functionally independent functions. 

A restricted version of this result was proven in ^ , where only the diagonal integrals la 
and the case A > was considered. However, the same algebraic results do hold for A < 0, 
and this possibility enable us to get other physical systems different from the one with A > 
that was solved in [Ij. We also remark that the existence of a (curved) Fradkin tensor Q is 
what makes T-Lx ^ a distinguished Hamiltonian, that is, a MS one which can be regarded as 
the "closest neighbour of nonconstant curvature" to the harmonic oscillator system, which is 
obtained in the limit A — )• 0. 

It is also worth stressing that theorem 1 can also be proven by relating Hx with the free 
Euclidean motion through a Stackel transform |21l [22| \2E[ [2^ [25] as follows. 

Let H be an "initial" Hamiltonian, Hu an "intermediate" one and H the "final" system 
given by 

H = ^ + V{cD, Hu = ^^+U{cD, H=^ = ^ + V{cD, (5) 
/x(q) ^(q) U ^(q) 

such that 

fl = fiU, V = V/U. (6) 

Then, each second-order integral of motion (symmetry) S of H leads to a new one S corre- 
sponding to H through an "intermediate" symmetry Sjj of Hjj. In particular, if S and Su 
are given by 

N 

S=Y1 ("i^P^Pj +W{ci) = So + Wiq) , Su = So + Wu (q) , (7) 
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then we get a second-order symmetry of H in the form 

Wrr 1 

S = So--^H + -H. (8) 

In our case, we consider as the initial Hamiltonian H ^ the free one on the NT) Euchdean 
space minus a real constant a (related with A and uj): 

H = ^p^-a, 2Xa = u}^. (9) 

And our aim is to perform a Stackel transform to the Hamiltonian T-Lx Q but written in 
"final" form as 



1 /p2 - 2a 

2 V 1 + Aq' 

Thus it can be checked that the transformation works provided that 



^ = ^A-a = -(i^^). (10) 



/, = 2, V = -a, ^ = 2(1 + Aq2), y = U = {l + W), (11) 

and the intermediate Hamiltonian is the A^D istropic harmonic oscillator 

^^C/ = ^P' + Aq' + I. (12) 

Next we consider the symmetries S of ([9]) which is clearly MS and endowed with 2 — 1 
functionally independent functions. Some of them are exactly ([3]): 

l<i<j<m 

•^M ^ X] ~ HPi)"^^ 'S'o,(m) = C(^rn), ^(m) = 0, 

N-m<i<j<N 

Sij=PiPj, So^ij=PiPj, Wij=0, (13) 



where m = 2, . . . , N and i, j = 1, . . . N . The symmetries Sjj of Hu (12) read 

Su,{m) = 'S'(m), 'S'o,(m) = C'(m)) ^l/,(m) = 0, 

'S'c/.ij = PiPj + 2Xqiqj, So^ij = piPj, Wu,ij = 2Xqiqj. (14) 



Consequently the Hamiltonian H (10) is also MS and its integrals of motion S ([8|) turn out 
to be 

Sim) ^ ^ (^.^ . _ ^.^^)2 ^ ^ ^ ^ 
l<i<j<m 

N-m<i<j<N 

Sij = PiPj - 2XqiqjH + H. (15) 

Finally by introducing H = T-L\ — a we recover all the results given in theorem 1 proving 
that, in fact, the Hamiltonian T-Lx is Stackel equivalent to the free Euclidean motion. 
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3 The underlying Darboux manifolds 



We recall that the real parameter A = 1/k was restricted in to take a positive value. 
Clearly, the superintegrability properties of the Hamiltonian stated in theorem 1 do hold for 
a negative A as well. Nevertheless the underlying space and the oscillator potential change 
dramatically with the sign of A in such a manner that the domain of the Hamiltonian must be 
restricted when A < 0. Hence, the "generic" Darboux space (that is, the Riemannian manifold 
with metric ^ determined by the kinetic part of ([l])) leads to three different manifolds A4^ 
which have the following geometric and topological properties. 



3.1 Type I: A > 

The Darboux space is the complete Riemannian manifold A4^ = (M^,g), with metric gij := 
(1 + Aq^) 5ij. The scalar curvature R{r) = ii(|q|) ^ is always a negative increasing function 
such that limr-!.oo R = and it has a minimum at the origin 

R{0) = -2XN{N-1), 

which is exactly the scalar curvature of the A^D hyperbolic space with negative constant 
sectional curvature equal to — 2A. 



3.2 Type II: A < restricted to the interior space 

In this case we consider the interior Darboux space defined by A4^ = {Br^,g) such that 

gij := (1 - |A|q2) 5ij, Br^ = [0,rc), r, = |q|c = 1/^, 

that is, Br^ denotes the ball centered at of radius rc which is the critical or singular value for 
which R{r) diverges and lim g,.^^^- = 0. It is clear that is incomplete as a Riemannian 
manifold. Notice also that 

R{0) = 2|A|iV(iV- 1), 

which coincides with the the scalar curvature of the A^D spherical space with positive constant 
sectional curvature equal to 2| A| . The behavior of R{r) depends on the dimension A*" as follows. 

• When 2 < A'^ < 6, the scalar curvature is a positive increasing function such that 
lim^.^^- R{r) = +oo. 

• If 7 < A^, there is a positive maximum for R{r) corresponding to 



N + 2 ^ _ 4\X\{N - l)iN - 2f 



y 2(Ar-2)|Ar -"-V-max;- 

and lim^^^- R{r) = -oo. 
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3.3 Type III: A < restricted to the exterior space 

To consider the exterior Darboux space, = {M.^\BT^,g), requires to change the sign of 
both the metric and scalar curvature Q: 

Qij := (|A|q2 - 1) 5ij, ^^XW, = (rc, oo), 
(iV- l)(2iV-3(Af-2)|A|q2) 



i?(q) = |A| 



(|A|q2-l)3 



Note that is again incomplete. According to the dimension A'^, the function R{r) behaves 
as follows: 

• For N = 2, this is a positive decreasing function such that lim^^^+ R{r) = +oo and 
lim,.^oo R^r) = 0. 

• If 3 < A < 5, the scalar curvature has a negative minimum 



A + 2 p. ^_ 4|A|(A-l)(A-2)3 



"^^^ Y2(A-2)|Ar ^ '""^ (A -6)2 

with lim^^^+ R{r) = +oo and limr_!>oo R{r) = 0. 

• When 6 < A, R{r) is a negative increasing function with lim^_j^^+ R{r) = — oo and 
limr._s.oo R{r) = 0. 

4 Three radial systems and their effective potentials 

Firstly, we remark that Tix can also be expressed in terms of hyperspherical coordinates r, 9j , 
and canonical momenta Pr,Pej, (j = 1, . . . , A — 1) defined by 

i-i N-i 
Qj = r cos 6j sin 

Ok, l<j<N, = r JJ sinOk, (16) 

k=l k=l 

SO, r = |q|. Thus the Hamiltonian ([T]) reduces to a ID radial system: 

= 1^ + ^(j^ = T(.,p„) (17) 



where \? = C^^^ = C(7v) is the total angular momentum given by 

iv-i i-i 

sm" 



iv-i j-i . 



i=l k=\ 



Now, the geometric analysis performed in the previous section indicates that we must 
deal with three different physical systems that, for the types I and II we name nonlinear 
hyperbolic oscillator and nonlinear spherical oscillator, respectively. In these two cases the 
generic expression for the Hamiltonian ([T]) is kept (with the boundary Tc for type II) , while for 



6 



type III the sign of the Hamiltonian has to be reversed, thus ensuring a positive kinetic term 
(and provided that the corresponding restriction on the domain is considered). In particular, 



as far as the nonhnear radial potential U{r) (17) is concerned we point out the following facts: 



• Nonlinear hyperbolic oscillator. When A > 0, the potential is a positive increasing function, 
such that 

m = ^^^y m = 0, and mnZY(r) = g. (19) 

• Nonlinear spherical oscillator. If A < and r < Tc, the potential is also a positive increasing 
function verifying 

^(r) = m ^(0) = 0, and lim U(r) = +oo. (20) 

2(l-|A|r^) r^rc- 

• Exterior potential. When A < and rc < r we impose the change of the sign of the 
Hamiltonian. In this way the potential becomes a positive decreasing function: 

^(^)= 9nxTy IV lim^(r) = +oo, lim Z^(r) = (21) 

2(|A|r^ — 1) r-!>rc+ r^oo 2\a\ 

But it is essential to stress that each of the above potentials has to be considered on 
the corresponding curved space described in section 3. In this respect, the complete classi- 
cal system can be better understood by introducing an effective potential (EP) that takes 
into account each curved background. This can be achieved by applying a ID canonical 
transformation |29| 

P = Pir,pr), Q = Qir), {Q,P} = 1, 



on the ID radial Hamiltonian (17) yielding 

nx{Q,P) = ^P^+U,s{Q). 
Next we present such an effective potential for the three abovementioned systems. 

4.1 The nonlinear hyperbolic oscillator 

The ID canonical transformation is defined by 

^(^■f')=7T^' «0 = ir^/TT^+=^^<^, (22) 
yl + Ar"' 2 2vA 

which implies that Q{r) has a unique (continuously differentiable) inverse r(Q), on the whole 
positive semiline, that is, both r,Q G [0, oo); note that dQ{r) = a/I + Ar^dr. This transfor- 
mation yields the EP 

Cat Up'r'^ 

where c^v > is the value of the integral of motion corresponding to the square of the total 



angular momentum C(Ar) = L (18). Hence the radial motion of the system can be described 



as the ID problem given by the potential Z^eff(Q('^))- 
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Hyperbolic/Spherical Uf.{{ 
35 r 




Figure 1: The effective nonlinear hyperbolic and spherical oscillator potentials ( |23| ) and ( |29| ) 
for A = ±0.02, Cat = 100 and uj = 1. The minimum of the hyperbolic potential (red curve) 
is located at = 3.49 with Wefr(^"mm) = 8.2 and Ues{oo) = 25, while the minimum of the 
spherical one (blue curve) is located at rmin = 2.86 with Z^eff ('^min) = 12.2 and rc = 7.07. The 
dashed line corresponds to the effective potential of the harmonic oscillator with A = with 
minimum Z^efr('^min) = 10 at rmin = 3.16. 



In fact, is always positive and it has a minimum located at rmin such that 



\CN + J X'^C% + Uj'^CN I 

- ^ , UM{r^ir.)) = -\CN + JX^C%+U^CN. (24) 



Therefore, rmin and Z//cff (Q(^'min)) are, respectively, greater and smaller than those corre- 
sponding to the isotropic harmonic oscillator, which are 



A = 



This EP has two representative limits: 



Cat/u;, ^efT(<5(?'min)) = ^^/cn- 



limZYeff (Q(r)) = +00, lim U^s{Q{r)) = uj^/{2\), 

1 — ^0 r— >oo 



the latter being coincident with (19). Thus, this EP is hydrogen-like (see fig. 1). 



(25) 
(26) 



4.2 The nonlinear spherical oscillator 

In this case, the canonical transformation is given by 

Pr 



TD( \ ni \ ^ R TTTT , arcsm(y A r) 

^(^.Pr) = , ^ , (5(r) = -rVl - |A|r^ + —= 

A/l-|A|r^ 2 2y^|A| 



so that Q(r) has a unique inverse r(Q) on the intervals 

1 



r E [0, rc 



TT 



(27) 



(28) 
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The EP reads 

which is always positive and it has a minimum located at rmin such that 



-| Alcjv + a/A^C^ + W^CAT 
Vin = J2 ' ^cff((3(?'min)) = |A|cAr + ^J A^C^ + Uj'^Cn- (30) 



But now rjnin and Z^cfr(Q(?'min)) are, respectively, smaller and greater than those correspond- 



ing to the isotropic harmonic oscillator (25). This EP has again two characteristic limits: 



limZYefr(Q(?')) = +00, lim ZYefr(<3(?')) = +00, (31) 

which means that we have a deformed oscilator potential that goes smoothly to an infinite 
barrier as r approaches rc~ (see fig. 1). 

4.3 The exterior potential 

The canonical transformation for the third system turns out to be 



1 , In 2 lAlr + v^VIAIr^-l 

P{r,pr) = , Q{r) = -rV|A|r2- 1 — ^ (32) 

and Q{r) has a unique inverse r{Q) on the intervals 

1 ln(2(v^|r; 
rG[rc,oo), Tc = Qg[Qc,oo), Qc = ^ ■ (33) 

The EP is ^ ^ 

z^eff(QW) = ^p^Sr^F + 2{\x\J^-iy 

The function lAeS is again positive but, unlike the two previous systems, it has no minimum; 



this fulfils the same limits (21 ) so EP is an infinite (left) potential barrier which is represented 
in fig. 2. 

Finally, some remarks concerning the quantization of these systems are in order. The 
nonlinear hyperbolic oscillator (A > 0) has been fully quantized in and its discrete 
spectrum is given by 



En = -h^\{n+^^ +h{n+^^^h^\^{n+^^ + (35) 

The corresponding stationary states have been obtained in analytic form. Note that the limit 
ra — )■ 00 of -En is just the asymptotic value u;^/2A, as expected. 

In view of the shape of the effective potential (see fig. 1), the quantum spherical os- 
cillator (A < 0) should provide a new radial confinement model that could be useful as a 
position-dependent-mass model for spherical quantum dots |30j . The exact solution of the 
corresponding Schrodinger problem is still in progress. 
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Figure 2: The effective nonlinear exterior oscillator potential (34) for A 
and u) = 1. The critical point is Vc = 7.07 and V(eg{oo) = 25 
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